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Abstract 

Relations between the gaugino masses have been shown to alleviate the degree 
of fine-tuning in the MSSM. In this paper we consider specific models of super- 
symmetry breaking with gravity mediation and demonstrate that within both GUT 
and string constructions it is possible to generate these relations in a natural way. 
We have numerically studied the degree of fine-tuning in these models, including 
one-loop corrections, and have found regions of parameter space that can satisfy all 
known collider constraints with fine-tunings less than 20%. We discuss some of the 
phenomenological features of these models within the regions of reduced fine-tuning. 

1 Introduction 

Low energy supersymmetry [1] (SUSY) is a well motivated proposal for physics beyond 
the standard model. In particular, the Minimal Supersymmetric Standard Model (MSSM) 
[2] predicts unification of the gauge couplings (to within a few percent) at a scale Mx ~ 
10^^ GeV and, together with R-parity conservation, provides a good dark matter candidate 
in the form of a neutralino. More significantly, perhaps, by stabilising the hierarchy between 
scales it explains why Mz <^ Mx and, through the soft breaking of supersymmetry, it may 
offer an explanation for the observed scale of electroweak symmetry breaking (EWSB). 
However, despite extensive searches performed at particle accelerators and dark matter 
detection experiments, no direct evidence for supersymmetry, or the Higgs, has yet been 
observed. 
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Following the absence of these signals at LEP, a lower bound of 1 14 Ge V was placed on 
the mass of the Standard Model (SM) Higgs, at 95% confidence [3]. Furthermore, limits 
have been set upon the masses of the supersymmetric particles, most notably the charginos, 
with0 > 88 — 102 GeV [1] . Several studies have shown that these constraints, combined 
with assumptions regarding the mediation scale and the pattern of soft terms, imply that 
a significant fine-tuning of the MSSM parameters is required if the correct scale of EWSB 
is to be obtained [3 El El E]. 

As is well known, it is the LEP limit on the SM Higgs mass that is the dominant source 
of this fine-tuning. Although the LEP limit does not strictly apply to the MSSM, it is 
relevant for regions of parameter space that contain a Higgs scalar with SM-like couplings. 
If CP- violating phases, 6i, are present in the soft terms, then this is not a generic feature 
of parameter space. However, the bounds on the electric dipole moments (EDMs) of the 
electron, neutron and muon place very strong constraints on these phases p]. Without 
some mechanism to suppress the EDMs, these phases are required to take small values and 
in the remainder of this paper we shall restrict ourselves to the case ~ or vr. In the 
absence of these phases, a SM-like Higgs is somewhat difficult to avoid. In the decoupling 
limit, where the CP-odd Higgs has a mass rriA ^ Mz, the lightest CP-even Higgs, h^, is 
SM-like. Its mass, M^o, including the dominant one- loop corrections from stop-top loops, 
is given by: 

Ml> = M|cos^2/3 + ^(in(^) + 5.) + ■ ■ ■ (1) 

which depends upon the top Yukawa coupling, hf, the top mass, Mt and the average stop 
mass M| = mi^mi^ (where m^-. are the eigenvalues of the stop mass matrix). M^o also 
depends upon a threshold correction, 5f, which is determined by the off-diagonal entry of 
the stop mass matrix Xt = At — /icot /3. From Eq ([1]) it follows that, for cos^ 2/3 ^ 1, the 
LEP limit requires a minimum stop mass: 

2 (T\/r2 _ t,j2\ 



27r2 (M^o - M, 



Ms > M^xp ^-, '^-^'\^r'' - 6t (2) 



Shf Mt 



which depends sensitively upon both ht and 6t. Previous numerical studies ^Q\, which 
include the two-loop 0{af, atas) corrections, have shown that for Mt = 173GeV and (5t ~ 
{St ~ 1) one typically requires Ms > 1 TeV (600 GeV). One can try to lower the bound on 
Ms by generating a large threshold correction, although this comes at the cost of requiring 
large values for At, and generally saturates at values 6t = Sf^^^ ~ 0{3). One can also 
consider the non-decoupling limit, in which rriA ~ Mz, where the two CP-even Higgs scalars 
are comparable in mass and involve large mixtures of the gauge eigenstates. Through 
this mixing the lightest Higgs can have a suppressed coupling to the Z, thus permitting 
M/jO ^114 GeV without conflict with LEP. The heavier Higgs, however, becomes SM-like 
in this case and, as pointed out by Dermisek and Gunion [11], heavy stops are still required 
to lift this Higgs above the LEP bound. Hence, in the absence of CP-violating phases at 
least, the LEP bound imposes significant constraints on the stop sector. 



^Where the precise limit depends upon the masses of the hghtest neutrahno, chargino and sleptons. 
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In models where the supersymmetry breaking is mediated at a high scale A ~ 0{Mx) 
these constraints have important implications for the scale of EWSB. Through radiative 
corrections the EWSB scale is related to the soft supersymmetry breaking terms [12] (our 
conventions for the soft supersymmetry breaking terms are described in Appendix A) and 
in general 

Xv^ = c^m2(A) + c'''Mj{A)Mk{A) + c'^V(A)-Mfc(A) + cV(A) (3) 

where v"^ = + whilst mf (A), A^j(A) and /i(A) are a scalar soft mass, dimension- 

1 soft term (either a gaugino mass, Mj, or trilinear term, Ai) and the supersymmetric Higgs 
mass, respectively, all defined at the mediation scale. The Higgs quartic coupling. A, is 
given at tree-level by A*-"-' = j{g'^ + gy)- The dominant radiative corrections are again 
generated by the top-stop sector, with contributions to f ^ of the form 

^ S 1^ (^^) K3(A) + mUA) + A?(A)) (4) 

Because of the high mediation scale A and ht ~ the soft terms associated with the 

stop sector appear in Eq ([3]) with 0{1) coefficients. For the Snowmass point SPS la [21], 
with A = Mx, one finds 

^ Xv^ ^0.3ml3{Mx) + 0.3mlj^{Mx) + l.5Ml{Mx) + 0.lA^t{Mx)-fi\Mx) + --- (5) 

The problem of fine-tuning arises because the constraint given in Eq ([2]) requires soft 
terms that generate contributions to Mz that are far too large. As a result cancellations 
must occur between different soft terms, and this cancellation requires fine-tuning of the 
parameters at the mediation scale. 

In order to measure the degree of fine-tuning we will follow previous studies of the 
MSSM [131 El E] and take the fine-tuning at a point in parameter space to be A = 
Max{|Aj|}, where Aj is the fine-tuning with respect to a parameter Oj, defined as: 

am ttj 

Values of A ~ 0{1) indicate that the EWSB scale arises naturally in that region of 
parameter space, with no special tuning of the parameters required. 

From the point of view that the MSSM is an effective theory below the mediation 
scale, each individual soft term is a parameter that should be included in the fine-tuning 
measure. Hence taking, for example, large scalar stop masses mq^^mu^ ~ 0{1 TeV) to 
achieve Ms ^ 1 TeV generates a fine-tuning Aq3 ~ 70 3> 1. This large fine-tuning 
indicates that a cancellation must occur between different parameters. In this effective 
theory picture, where all soft masses are unrelated, such a cancellation would be unnatural. 
The probability p of such a cancellation is often taken as p = A^^ ~ 1%, which implies 
that the observed EWSB scale is a product of chance in this region of parameter space. 

This fine-tuning may, however, be an artifact of the effective theory description. In 
a more complete theory of supersymmetry breaking the soft terms of the low energy La- 
grangian will be related and, written in terms of the parameters of the true theory, the 
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observed EWSB scale may prove to be natural with no fine-tuning at all. Previous studies 
of the MSSM have tried to identify relations that, if present amongst soft terms at the 
mediation scale, significantly reduce the level of fine-tuning (TJ El [12]. Most interestingly, 
some of these relations can be easily motivated within known models of supersymmetry 
breaking. 

One such relationship is the universality of the scalar masses at the high scale: 

m^{Mx)=ml (7) 

which is a natural consequence of a flavour-blind mediation mechanism, such as supergrav- 
ity. In this case the EWSB condition of Eq ([3]) becomes: 

Xv^ = d^ml + d^Mi{Mx)Mj{Mx) + ■■■ (8) 

It has been shown that, provided tan /3 > 0(5), a cancellation occurs between the tree 
level soft masses and the radiative corrections, such that c™ <^ 1 [Tll|15]. This reduces the 
dependence of Mz upon the parameter mo, hence permitting A^^ ~ 2c'"mo/M| < 10 even 
for large mo ~ 1 TeV. This permits the heavy stops required to satisfy the LEP bound on 
the Higgs without a large fine-tuning. 

Relations between gaugino masses have also been shown to reduce fine-tuning. In 
some of the earhest fine-tuning studies [8] it was noted that Mz is strongly dependent 
upon the gluino mass M^{K), thanks to the large values of both ht and as- Hence it 
was suggested that a small fine-tuning lAi/gl < 10 requires a relatively light gluino, with 
M^{h) < 160 GeV, whilst the bino and wino masses may be much heavier. Later studies 
have since pointed out that relations between gaugino masses can permit a much heavier 
gluino without significant fine-tuning [16]. The basic idea is that the gaugino masses are 
related to some common parameter, such that Mj(A) = riiMi/2, where the rji are fixeci 
ratios. In this case the contribution of the gaugino mass parameter M1/2 to the EWSB 
scale is given by: 

Xv^ = c'm^iA) + c^'Ml/^ + ■■■ (9) 

and for some choice of the rji the co-efficient c^^ <^ 1. Hence, one can have ~ -'-0 

even with M1/2 ~ 0{1 TeV). This idea was discussed recently by Abe, Kobayashi and 
Omura (AKO) [19], who identified, in a bottom-up approach, certain non- universal ratios 
which permit heavy gaugino masses without a large fine-tuning. Furthermore, they point 
out that these large gaugino masses can, through radiative corrections, generate the heavy 
stops needed to lift the Higgs above the LEP limit. 

Of course the critical question is whether the gaugino mass ratios that reduce fine- 
tuning can be generated by a specific model of supersymmetry breaking, otherwise the 
initial choice of gaugino mass ratios is itself fine-tuned. In this paper we consider in detail 
specific supersymmetry breaking models with gravity mediation that do reduce fine-tuning. 
The layout of our paper is as follows. In Section[2]we will identify a preferred set of gaugino 
mass ratios that will reduce fine-tuning. Some of this work will overlap with the study of 

^By fixed, we mean that the rji are not continuous parameters. Instead they are discrete ratios, generated 
by some mechanism, and because of this they are not included in the fine-tuning measure. 
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AKO, but we include it here for completeness. In Section [3] we will compare these ratios 
to those that are predicted by known mechanisms of generating non-universal gaugino 
masses in supergravity. For those models that predict ratios in the preferred range we have 
performed a detailed numerical study of the fine-tuning, which we present in Section |H In 
Section |5] we highlight some of the phenomenology of these models in regions of parameter 
space with A < 10. Finally, we conclude in Section |6l 



2 Reduced fine-tuning through non-universal gaugino 
masses 

In this section we identify the relations between gaugino masses that reduce fine-tuning. 
We will work within the context of the R-parity conserving MSSM, and we will assume 
no additional field content below the unification scale Mx- We will assume that the 
supersymmetry breaking is mediated to the observable sector at a high scale A. In general 
A, and the field content between the scales Mx and A will depend upon the specific 
details of the model under consideration. In this paper we will make the simplifying 
assumption that A = Mx, and we will ignore the effects of threshold corrections due to 
states lying above the scale Mx- Furthermore we will assume that some mechanism, which 
we will discuss further in Section El generates gaugino masses that are related to a common 
parameter M1/2: 

M,{Mx) = n-Mi/2 (10) 

where the ratios rji are either fixed or discrete, i.e. cannot be continuously varied. Hence 
whilst the parameter M1/2 must be included in the fine-tuning measure, the r]i are not. 
Without loss of generality we can define M2{Mx) = M1/2 (i.e 772 = 1) and we can always 
choose our phases such that M2{Mx) is positive. We will also assume that the scalar 
masses and trilinear terms adopt universal relations: 

m," iMx)=ml Ai,jk{Mx) = Ao (11) 

Hence the continuous, dimension-1 parameters that must be included in the fine-tuning 
measure is the set a = {mo, M1/2, Aq, B{Mx), fi{Mx)}- We will refer to a generic member 
of this set using a subscript, a^. 

In order to calculate the fine-tuning, it is necessary to determine the relationship be- 
tween V and the parameters a,. After minimising the one-loop effective potential one 
obtains a relation between v and the running soft terms of the form: 

2 ^ tan2^ _1 _ 2 ^^2) 

tan2/3-l tan2/3-l ^'^ ^ ^ 

where A*^"-* = j{g'^ + gy) is the tree-level quartic coupling and rn^^ = m'jj^ + dV^^^ dH^. 
V^^^ is the one-loop Coleman- Weinberg (CW) potential pjj. As usual, we define tan/3 = 
{H^} I (H^) , which is also related to the soft terms: 



sin 2/3 = / ^, „ . 13 
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The dominant 0{hf) contribution to the CW potential is generated by the top-stop 
sector: 



dm 



3hl 



Q3 



m^3 + X^) In 



[m 



ti 



m 



t2' 



t2' 



In 



t2 



_2/iyin(^*^ 



(14) 



where Q is the renormahzation scale and the m^- are the eigenvalues of the stop mass ma- 
trix. Note that all of the quantities appearing in Eqs (!T4|) and ( [T2|) are running quantities, 
with an implicit dependence upon Q. In order to correctly resum large logarithms and 
minimise the corrections from the CW potential, it is evident from Eq ( !T4l) that one should 
choose a scale Q ~ Ms = ^fnifjn^. In the remainder of this paper we will fix Q = M5, 
unless stated otherwise. By solving the RGEs [20] the running soft terms at the scale Ms 
can be related to the supersymmetry breaking parameters ai. 



m^Ms) 

Ma{Ms) 



zi'^ajak 



MA ; 



'M 



(15) 



B[Ms) = B{Mx) + Z^'My2 + Z^A, 



whilst the /x parameter is given by: 



^x{Ms) = Z^fxiMx] 



(16) 



The coefficients zj^ and have an implicit dependence upon the scale Q = Ms and the 
boundary conditions at Mx- Using the expressions for the running quantities given in 
Eqs f lTSj) and ( !T6|) one can then express Eq (fT2|) in terms of the parameters that describe 
the supersymmetry breaking. This gives a relation between these parameters and v. 



/2 



c^'^My^Ao + c^Al + c^/i^ + c*'^Mi/2/x + c^Mo/i (17) 



where A = A*^°^ + A*^^\ and A*^^) is a one-loop contribution to the quartic coupling. It should 
be clear that the coefficients c and A, although dimensionless, depend upon the parameters 
Qi and V through the structure of the one-loop corrections and through the dependence of 
tan^ upon the soft terms, which is given by Eq f|T3l) . 

Let us now consider the fine-tuning with respect to a parameter a^. For the measure 
given in Eq ([6]) this is: 



al dMl 



2A(°) 



where in the approximate equality we have neglected the self-energy of the Z boson, 11^^, 
that contributes to its mass through M| = 2\^^'^v^ + Re Yizz- Using Eq ^ it is straight- 
forward to calculate the derivative dv"^ /daf and determine the fine-tuning. However, let us 
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note that for the choice of renormahzation scale Q = Ms the contributions from the CW 
potential are suppressed by loop factors of (47r)^ with no large logarithmic enhancement. 
Hence, to begin with, let us ignore these one-loop corrections and calculate the fine-tuning 
at the leading tree level. Furthermore, since we want to minimise the degree of fine-tuning 
we will restrict ourselves to the region of parameter space in which tan^ /3 ^ 1. As can be 
seen from Eq ([1]) this maximises the tree level contribution to the Higgs mass, and thus 
alleviates the constraints on the stop sector that result in the fine-tuning problem. In this 
case it follows from Eq f[T7|) 



= -z]}ml - zlM!/, - zlAl - 2zl^My,Ao - Zy{A) + ■ ■ ■ (19) 

where the ellipses imply terms that are suppressed by factors of tan^ (3 or one-loop coeffi- 
cients (47r)^. The leading contribution to the fine-tuning A^, due to the gaugino mass, is 
therefore given by: 

zfl Ml, + zff^Mu^Ao 
Am ^ -2 ^" J" (20) 



z 

It is clear from Eq ( 120|) that if the coefficients z^^ and z^^ are suppressed, then one 
can have a small fine-tuning | Aj\/| < 10 even with a large gaugino mass Mi/2- For example, 
one can achieve \!\m\ < 10 for M1/2 = 1 TeV provided zjj[ < 0.04 ~ (87r)-\ if Aq = 0. 
The ratios rji that therefore reduce fine-tuning are those for which z^^ is suppressed. 

In order to calculate the RG coefficients of Eq f|T5l) and their dependence upon the rji it 
is necessary to specify the boundary conditions at the unification scale. In particular, we 
need to know the value of Mx together with the values of the gauge and Yukawa couplings 
at that scale. These couplings must be determined by matching them onto low energy 
data, including threshold corrections from superpartners. To perform this calculation 
we have used the software package SOFTSUSY 2.0 [21J, together with the following low 
energy inputs: asiMz)'^ = 0.1176 ± 0.002, a-\Mz)'^ = 127.918, Mz = 91.1876 GeV, 
rribimh)'^ = 4.20 GeV [22] and Mt = 172.6 ± 1.4 GeV [2S]. Shown in Tableware the high 
scale couplings obtained for a number of Snowmass benchmark points [24] . Note that the 
errors in Mx, ht{Mx) and g^^Mx) are the propagated errors due to the la uncertainties 
in the top pole mass Mt and as{Mz), added in quadrature. 

Because the couplings depend upon threshold corrections they are functions of the 
spectrum of superpartners, and hence of the parameters Oj and rji. However, we note from 
Table [Tjthat the uncertainty due to these unknown threshold corrections is relatively small; 
comparable to the uncertainty due to errors in the low energy data. Thus, in what follows, 
we will fix these high scale couplings to the values appropriate for one supersymmetric 
spectrum, and then neglect further threshold dependence as we adjust the ratios rji and 
the parameters a^. 

Using the central values of SPS point la we find that the coefficients zff^ and zff^ have 
the following form at a scale Q = 500 GeV: 



zfj[ = -1.53r/| + 0.21 - 0A3r]3 - 0.020r/ir/3 + 0m7r]l - O.OOSt^i (21) 

MA 



z;^^ = 0.29r/3 + 0.07 + O.Olr^i (22) 
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Benchmark Point 


Mx 1 10^6 


htiMx) 




gAMx) 


la 


2.39(2) 


0.494(18) 


0.708(2) 


0.721 


2 


2.38(2) 


0.506(18) 


0.703(2) 


0.714 


3 


1.97(2) 


0.498(18) 


0.705(2) 


0.716 


5 


2.15(2) 


0.508(20) 


0.707(2) 


0.718 



Table 1: Values for Mx and the couplings calculated using SOFTSUSY for four different Snow- 
mass points. The uncertainty is calculated by allowing Mt and as{Mz) vary over their one a 
range, and then combining the resulting uncertainties in quadrature. The uncertainty in gi is 
small, typically less then 0.1%. 



It is clear that the contributions from the gluino mass are the most significant, as they are 
generated by the larger, and faster running, couplings ht and a^. Note also that, in zf^^, 
the coefficient of r^f is negative because the gluino mass feeds into the evolution of mfj^ 
through a Yukawa coupling. This is critically important, since it permits a set of ratios 
rii such that zff^{Q = 500 GeV) = 0. We note that there is no set of ratios such that 
both zff^ = zfl^ = 0. For |?7i| < C(l) we observe that the coefficient zff^{Q = 500 GeV) 
vanishes for 773 ^ 0.33, —0.42 , whilst z^^ ^ 0.17, —0.05 respectively. 

From this discussion it should be clear that for any choice of ratios rji the coefficient 
zff^{Q) will vanish at some scale. At this scale the running soft mass m^Hu ^^^^ become 
independent of the parameter M1/2, provided Aq vanishes. This is analogous to the scalar 
focus point scale [12], where m'jj^ becomes independent of the scalar mass parameter mg. 
From the perspective of an RG analysis it is this scalar focus point, and its proximity to 
the TeV scale, which is responsible for the smallness of the coefficient appearing in 
Eq ([8]). Because of this analogy we will call the scale at which zff^ = the gaugino focus 
pointo, Q^^ . In Fig [1] we plot this focus point scale in the i]^ — r/i plane. 

For Ms ~ Q^'^ the coefficient z^^{Q = Ms), that appears in the fine-tuning Am, has 
the following form: 



Hjf) (23) 

-.QM Ms 



where the coefficients zf^ and that appear in Eq fl2^ are given in Appendix B. Thus, 
provided is sufficiently close to Ms the coefficient z^^ can be suppressed so as to permit 
large values of M1/2 without a significant fine-tuning. Because the coefficient z^^ cannot 
be simultaneously small, |Am| < 10 will require that Aq is small for large Mi/2- Hence we 
will restrict ourselves to the plane = for the remainder of this papei0. 

As we discussed in Section [T] the LEP bound on the Higgs places a lower bound on 
the stop mass Ms through Eq ([2]). For mg = = 0, one can determine the value of 



■^This ignores the minor technical point that the gaugino focus point involves cancellations between 
radiative corrections, whilst the scalar focus point involves cancellations between tree level and radiative 
pieces. 

''Abe et. al. have explored the possibility of allowing a large Aq by postulating a further relation of 
the form Aq = rj^Mx/^- 
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Figure 1: A contour plot of the gaugino focus point scale Q (GeV), as a function of the gaugino 
mass ratios r/i and 173. The light grey (dark grey, hatched) region indicates the ratios that will 
permit Am < 10 for M5 = 600 GeV (M5 = 1 TeV). The filled circles, squares and triangles 
indicate the mass ratios predicted by the SU{5), 5*0(10) and Eq GUT models, respectively, that 
are discussed in Section [3l The stars give the prediction of a string model with moduli-dominated 
SUSY breaking (the O-II model |29]), for various integer values of the parameter 5gs- The ratios 
that are permitted within models with mixed Moduli- Anomaly breaking lie upon the dashed line. 



Ml/2 = M which achieves this minimum stop mass for a given set of ratios rji. This 
is straightforward to calculate using the RG coefficients given in Appendix B. Moreover, 
|A^| = |1 — Am I < I Am I in this tree- level approximation, thus A = |Am|- The fine-tuning 
introduced by the LEP bound is then, for a particular set of ratios rji, 



A 



-H. 



[{Ms)M 



Ml 



(24) 



In Fig [T] we have highlighted, in the shaded and hatched regions, the ratios that have 
A < 10 for Ms = 600 GeV and Ms = 1 TeV, respectively. These values of Ms, according 
to a previous numerical study [lOj, are sufficient to lift the Higgs above the LEP bound 
when 5j ~ 1 and 0, respectively. The ratios in these shaded regions are therefore the 
'preferred' ratios that reduce fine-tuning. An interesting question to ask is whether any 
known models predict such ratios, and we will turn to this in Section |3J 

It is worth making a number of comments about this analysis. We have investigated 
the effect of uncertainties in the low energy data, namely in Mt and as{Mz), combined 
with uncertainties associated with the unknown threshold corrections. Allowing the low 
energy variables to shift by their la errors and also choosing different Snowmass benchmark 
points, we observe only minor shifts (< 1%) in the position of the preferred region. What 
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is more important is the effect of calculating A beyond the tree level approximation. For 
the ratios with a focus point ~ Ms it is clear that will be comparable to one-loop 
contributions. Therefore, in the preferred region of FiglH one expects that the fine-tuning 
given by Eq f lM|) is only an order of magnitude estimate. Thus the results presented in 
Fig [1] should be used only as a guide to which ratios deserve a more detailed study. To 
determine the fine-tuning more accurately it will be necessary to calculate all of the one- 
loop and tan /3 suppressed contributions. Furthermore, the RGBs should be solved at the 
two-loop level, since the sub- leading logarithms arising at two loops are comparable to the 
one-loop contributions from the CW potential. We will perform such a study, numerically, 
for any model that predicts ratios in the preferred region. 



3 Non-universal gaugino masses from GUTs and strings 

In the previous section we have identified a preferred set of gaugino mass ratios that, 
potentially, could alleviate fine-tuning within the MSSM. Hence it is clearly of interest 
to examine whether any known models of supersymmetry breaking predict these ratios. 
For simplicity, we will restrict ourselves to models where the supersymmetry breaking is 
mediated by gravity alone. In this case there are several known mechanisms for generating 
non-universal gaugino masses, which we will briefly review here. 

In general, the gaugino masses for a supergravity model are given by [25] : 

M^p = :^^e-^/'G\G'y,{df:,/d<f>n (25) 

where G is the Kahler potential, a real gauge singlet function of the chiral superfields 
and fap is the gauge kinetic function, an analytic function of the $i that transforms under 
the gauge group as the symmetric product of two adjoints, {Adj x Adj)s- From ( 125|1 it 
is clear that non-universal gaugino masses will result if fap has some non-trivial gauge 
structure. 

In a GUT this structure can be generated by a Higgs field, E, provided it transforms 
under the gauge group G in an irreducible representation R C {Adj x Adj)s- If this field 
takes a vacuum expectation value (VEV), spontaneously breaking G into a subgrouplf] H, 
it can produce a gauge non-singlet contribution to /q,^. In particular, for a gauge kinetic 
function of the form: 

U = S^p + B{<t>,)^ (26) 

where B is an arbitrary (gauge singlet) function of the $i, the gaugino mass matrix is 
given by [26] : 

Map = vJapM{T.) (27) 
^Where this subgroup must include the Standard Model. 
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G 


R 


H 


V3 


Vi 


SU{h) 


75 


SU{3) X SU{2) X U{1) 


1/3 


-5/3 




54 


SU{2) X S0{7) 


-3/7 


-3/7 


50(10) 


210 


[SUibY X t/(l)]flipped 


1/3 


-1/15 




770 


[su{5y X f/(i)]flippcd 


-2/3 


101/15 



Table 2: The ratios of gaugino masses generated at Mx for the subset of SU{5), SO{10) and Eq 
GUTs [261 EH EH] that appear in the preferred, shaded region of Fig [TJ 

where the discrete coefficients rja are determined by R and H. For the GUT groups SU (5), 
50(10) and Eq the possible ratios have been determined and can be found in the literaturqj 
[261 EH [28]. 

In Fig d] we have plotted the focus point scale for the ratios rja predicted by the GUT 
models. From the figure we note that there are several models that have in the region 
that may permit low fine-tuning. The relevant details of these models are given in Table [21 
and we will study these in more detail in Section [H As pointed out by Martin [28], the 
54 breaks 5*0(10) to a subgroup, SU{2) x 50(7), which cannot be reconciled with the 
chiral fermion content of the SM. Hence the 54 is not a realistic model. However, it does 
generate mass ratios that illustrate an interesting phenomenological structure and, pending 
the identification of a viable model with mass ratios in this region, we will include it in our 
phenomenological survey. 

String theories also provide a mechanism for generating non-trivial structure in fajs, in 
the form of contributions related to anomalies and string threshold corrections [29l [30] . Al- 
though these arise as one loop corrections to the gauge kinetic function, they can generate 
large non-universality in the gaugino masses provided that the supersymmetry breaking 
occurs in the moduli, and not the dilaton, fields. The form of the string threshold correc- 
tions is highly model dependent, as they are sensitive to a number of factors, including: the 
choice of manifold for the compactified dimensions, the modular weights of the observable 
sector fields and the moduli VEVs. It is well beyond the scope of this work to survey all 
the possibilities; here we will only consider a simplified model based on an orbifold com- 
pactification. This model, coined the 'O-H' model [29], generates ratio^r^Q oc {—Scs + ba), 
where ba = (33/5, 1, —3) is the one-loop beta function coefficient for the gauge couplings, 
whilst 5cs is a constancl that is required to be a negative integer by an anomaly cancella- 
tion condition. Since 5gs is a discrete parameter it need not be included in our measure 
of fine-tuning. 

Non- universal gaugino masses can also be generated in the so-called 'Mixed Moduli- 

^The recent study by Martin [28] has highlighted possible errors in the previous study of S'O(IO) [27] , 
In this work we use the results of Martin. 

''Note that, unlike the case explicitly considered in [29], we assume VEVs for the moduli fields that 
are far from the self-dual point. In this case the contribution from the superconformal anomaly can be 
neglected [30] and a large hierarchy between the gaugino and scalar masses can be avoided. 

®This constant is associated with the Green-Schwarz counter term that must be introduced to cancel 
anomalies. 
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Anomaly' (MMA) scenario [31j, where the moduh fields and superconformal anomaly gen- 
erate comparable contributions to the soft terms. In such models the gaugino mass ratios 
have the same structure as in the 'O-II' model: rja = a + b^, where a is the ratio of the 
moduli and anomaly contributions. In this case the parameter a is related to topological 
quantities, and is a rational number involving the ratios of fiuxes. Despite the fact a takes 
only discrete values, if the interval between these values are small, it should be included in 
the fine-tuning measure. In this case there is no reduction in the fine-tuning. This should 
be compared with the O-II model where the Green-Schwarz term takes on only integer 
values and thus does not contribute to the fine-tuning measure. 

The mass ratios predicted by the MMA scenario lie along the dashed line shown in 
Fig [H whilst those of the O-II model are indicated by stars. We note that the O-II 
with 6gs = ~5, which has gaugino mass ratios {rji = 29/15, 773 = 1/3}, possesses a set of 
ratios that may permit low fine-tuning. In fact, this point produces the same gaugino mass 
ratios predicted by the 'Mirage Mediation' point of the MMA scenario, which has previously 
been demonstrated to possess low fine-tuning [31j. This set of ratios is characterised by the 
unification of the gaugino masses at a low scale. Because the MMA scenario is complicated 
by the need to assign modular weights to each of the fields, we will postpone its discussion 
to a later paper. We simply note here that a low scale gaugino focus point can, for an 
appropriate choice of a, be achieved in the MMA with M3M2 < 0. 

In order to make more quantitative statements about the naturalness of these models 
we must calculate the fine-tuning beyond the tree level approximation. We have done this 
numerically and we discuss our results in the following section. Note that, henceforth, we 
will refer to the O-II model with 6gs = —5 as simply the 'O-II', and we will refer to the 
GUT models by the representation R that is listed in Table [H 

4 Fine-tuning in the O-II and GUT models 

Using SOFTSUSY, together with the low-energy inputs discussed in Section [21 we have 
calculated the fine-tuning in a fixed grid scan over the parameters mo and Mi/2- The 
derivative dMz/dai is calculated by SOFTSUSY using the following procedure: 

• A parameter a, is shifted by a small quantity at the scale Mx- 

• The couplings and soft-masses are run down to the low scale Ms using two loop 
RGBs. 

• The one-loop potential is then minimized to obtain the new Higgs VEV, v. 

• This procedure is then repeated, calculating v for increasingly smaller shifts of Oj. 

By default, this procedure neglects the v dependence of masses (of the fermions and 
scalars), appearing in the one-loop potential. In order to calculate the fine-tuning cor- 
rectly at the one-loop level we modified this procedure, so that the one-loop potential is 
correctly minimised including the dependence of all masses upon v. 
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We show in Fig |2] the results of this scan for the O-II, 54, 770 and 210, for Aq = 0, 
tan/3 = 10 and fi > 0. The case /i < is shown in Fig [31 Each plot contains several 
shaded regions: a hatched (red) region where the conditions for successful EWSB are 
not satisfied; a dark, unhatched (blue) region where the LEP constraints on the masses 
of superpartners are violatecl^ and a lighter, unhatched (green) region where the lightest 
stable particle is charged, which is in conflict with cosmological observations and searches 
for anomalously heavy molecules [33] . (We will discuss these last two regions in more detail 
in Section |3) The unbroken lines (black) are contours of constant fine-tuning. A, whilst 
the dotted lines (black) are contours of constant A^, which we have included so as to 'guide 
the eye', and assist in identifying the regions where A = |A^|. The dashed and dot-dashed 
lines (both shown in red) are contours of constant Higgs mass, with M^o = 111 GeV and 
M/jO = 114 GeV, respectively. 

As can be seen in Figs H] and [31 EWSB fails to occur in a significant region of the 
mo — Ml/2 plane. This is straightforward to understand by considering the structure of the 
radiative corrections. At the minimum of the one-loop potential the following condition 
must be satisfied: 



\(0) 2 , I |2_ m 2 , ^Ha^O + , 'irq [, 2 , 2 , 42m 

X^^'h ^-z mo + + I (mQ3 + m^3 + A J In 




(m?-mn + f4 f^Al -2M2ln,— 2 




where here we have only included the dominant 0(/i^) contribution. This condition is 
used to determine \fi\ at each point in parameter space. Note that we have chosen our 
renormalization scale Q = Q^^ such that z^' = 0, so that the dominant dependence of /j, 
upon the gaugino mass, M1/2, is shifted to the one- loop corrections. For this choice we 
find that z"^ is small and positive (~ 0.03), as the scalar focus point scale lies below Q^^ in 
each of these models. It is evident from Eq ( |28l) that as the stop mass scale M| = mi^mi^ 
increases and approaches eQ^^ the contribution of the one-loop corrections to Eq f[28|) will 
change sign. This will cause \ix\ to decrease with increasing M5, until eventually < 0. 
At this point EWSB can no longer occur. The 75 model, which according to Fig [H should 
also exhibit low fine-tuning, is a more extreme example. From Fig [Hit can be seen that this 
model has a focus point ~ 400 GeV, which is significantly lower than the other models 
considered here. Because of this, the 75 permits EWSB in only a very small region of the 
mo — Ml/2 plane. One can increase the size of this region by permitting non-zero values 
for Aq, which can drive to more negative values. We, however, will not discuss this 
possibility further here; we simply note that a low scale focus point is not always sufficient 
to guarantee satisfactory EWSB combined with low fine-tuning. 

Within the region that permits EWSB each model exhibits a fine-tuning. A, which is 
less than 10. Furthermore, there is a substantial region in which A < 5. This is to be 
compared with the CMSSM (which assumes a universal gaugino mass, rj^ = 1), which is 



shown in Fig 4(a) for tan /3 = 10 and Aq = 0. In the CMSSM the region with A < 10 



^In most cases it was the limit on chargino masses that proved to be the most constraining. 
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mo/GeV mo/GeV 
(c) 770 (d) O-II 

Figure 2: The solid (black) contours indicate the fine-tuning in the mo — plane. This is 
for the hypersurface in parameter space with tan P = 10, Ao = and > 0. The (red) dashed 
and (red) dot-dashed contours indicate where the Higgs mass, m^o, is 111 GeV and 114 GeV, 
respectively. The dotted contours are lines of constant A„. 



(A < 5) is mostly (completely) excluded by constraints from chargino searches at LEP. As 
discussed in Section |2l it is the low scale focus points for both mo and M1/2, present in 
the 54, 210, 770 and O-II, which suppresses the dependence of Mz upon these parameters 
and permits |Amol, lA^i/al ~ 10 these models even for mo, M1/2 ~ 1 TeV. However, the 
fine-tunings observed in Figs |2] and [3] are somewhat smaller than what one would naively 
expect from a tree level analysis. For example, at tree level the derivative dv'^/dfi'^ is given 

^^'^di:^^~'^^'i^H^ ^^^^ 
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Figure 3: These are the same plots as in Fig [2l but for the case < 0. The (blue) dashed 
contours, with long and short dashing, indicate the sensitivity to the top Yukawa, A/^^. The 
contours for A/^^ are labelled in italic. 

Hence the fine-tuning in the /i parameter is given by: 

which, neglecting the second term appearing in Eq ( l30ll . gives an upper bound of |yu| < 
200 GeV < 140 GeV) for regions with A < 10 (A < 5). Whilst we find that each 
model satisfies |yu| < 200 GeV over the entire mo — M1/2 plane, we find that > 140 GeV 
in regions with A < 5. There are two reasons why this is so. First, for the case of moderate 
tan/3 = 10, the second term0 appearing in Eq f l30|) is comparable to the first for a Higgs 

^°The presence of this second term is the direct result of assuming that the soft mass that mixes the 
Higgs fields is of the form BjiHuHd- 
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soft mass ~ 1 TeV. Hence a mild cancellation arises, which tends to reduce |A^|. 
Second, when the one-loop corrections are taken into account, the derivative dv^ ld\j? is 
given by: 

where e is given by: 

^^AU (pg^ , ,n C'l^) (2 - <pL±^ 1 .... 1 (32, 



and is related to the one-loop contribution to the Higgs quartic coupling. Hence the fine- 
tuning is reduced by an overall factor of (1 -|- e) ~ M^o/M^. In regions where the Higgs 
mass M/jO ~ 114 GeV, this reduction is a factor between 1.5 — 2. Thus the large stop 
masses required to lift the Higgs mass above the LEP bound assist in ameliorating the 
fine-tuning. 

The structure of the radiative corrections, combined with a low focus point scale, reduce 
the fine-tuning for the other soft parameters in another way. Consider the right hand side of 
Eq f l28|) as a function of the parameters aj and v (where here aj includes soft parameters 
only), 

A(°)t;2 + ^ J^(a„/i,t;2) (33) 

In the limit that the soft terms vanish and supersymmetry is restored, J-' must vanish. 
However, as we have noted above, J-" will change sign for sufficiently large Ms, because 
of the low scale focus points. Thus J-" has another zero for a non-zero value of the soft 
masses. Between these two zeros the function J-" will possess a maximum, at which point 
the derivative dv^/dai will vanish: 



and hence the fine-tuning will vanish also. We demonstrate this explicitly in Fig 4(b) which 
shows a plot of I Am I in the O-II model, for a scan over M1/2 with fixed mo = 500 GeV. 
As can be seen, for M1/2 ~ 1 TeV the fine-tuning vanishes. The tendency of the radiative 
corrections to reduce the level of fine-tuning is not a new observation, as it has been noted 
in several previous works. However, in the scenarios considered here, with low scale focus 
points in both the gaugino and scalar soft masses, the one-loop and tree level contributions 
to the fine-tuning are comparable, causing this effect to become pronounced. 

Let us now discuss the effect of the LEP bound on the Higgs mass. The theoretical 
calculation of the Higgs mass currently has an uncertainty at the level of 3 GeV[32]. As 
a result, we will only consider points in parameter space with M/jO < 111 GeV to be in 
violation of the LEP limit. It can be seen in Figs [2] and Elthat all of the models contain 
regions of the mo — M1/2 plane that can satisfy the LEP constraint, with fine-tunings A < 5. 
It is also clear that the LEP bound is far more constraining in the models with M3 < 0, 
that is the 54 and 770. This can be understood if one considers the size of the threshold 
correction 6t that, as shown in Eq ([T]), appears in the correction to the Higgs mass. This 
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Figure 4: Shown in Fig 4(a) is a plot of the fine-tuning contours in the CMSSM, for tan/3 = 10, 



= and ^ > 0. The shaded regions and contours are the same as in Fig[3l Shown in Fig |4(b)] 
is a plot of I Ajv/I as a function of M^j^^ for the O-II model with tan /3 = 10, Ao=0, mo = 500 GeV 
and fi > 0. 



threshold correction depends upon the degree of stop mixing and is determined by the off 
diagonal entry in the stop matrix Xt = At — fi cot /3. In the limit that the stop soft masses 
are equal, m^g = mfj^, 6t has the following simple form: 

5. = X,^(l-§) (35) 

where Xt = Xt/Ms- (Although the stop soft masses are not equal in the models we 
consider, this limit serves to neatly illustrate our point.) Because /icot /3 is typically small 
in these models (< 20 GeV), Xt ~ At. For the case At{Mx) = that we consider here, 
the trilinear term at the low scale is generated solely through radiative corrections. At a 
scale Q = 500 GeV the trilinear term is given by: 

At(Q = 500 GeV) = My2 ( - l-6l7?3 - 0.25 - O.OSr^i) (36) 

Hence in models with r^a < a cancellation occurs between the first and second terms, 
which tends to suppress the ratio Xt- As a result the 54 and 770 models have smaller 
values for 6t (we find that 6t ~ 0.2 and 0.1 for the 54 and 770 respectively), compared to 
the 210 and 0-11 (which both have 6t ~ 1.1). From Eq ([2]) it thus follows that they need 
heavier stops in order to satisfy the Higgs bound. 

Thus we have identified several models that can satisfy the LEP bound and give the 
correct EWSB scale without significant fine-tuning in the dimensionful parameters. How- 
ever, one may worry that the models contain other sources of fine-tuning. Here we will 
briefly discuss two possibilities: tunings in the dimensionless parameters and tunings that 
may arise at larger values of tan/3. 
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4.1 Tunings in ht{Mx) 

As pointed out by Romanino and Strumia [51], it is difficult to make a case for a model 
that claims to reduce fine-tuning in the soft terms if it requires a large tuning in its 
dimensionless parameters — in particular the top Yukawa coupling at the high scale ht{Mx)- 
This parameter is unlike the other dimension- 1 parameters, since its value is fixed by low 
energy observables and the spectrum of superpartners. For the SPS point la we find 
ht{Mx) = 0.494 ± 0.18, where the error combines uncertainties from the measured values 
of Mt and as{Mz)- To reflect this difference we follow previous works by calculating the 
fine-tuning in ht using the following measure: 

_5h, d\nM, 

- 17 m^TT^ ^^^^ 

where 6ht is the uncertainty in ht{Mx)- In Fig[3]we plot contours of A/^^, which are shown 
with (blue) alternating long and short dashes, assuming that 6ht/ht = 4%. This fractional 
uncertainty is slightly larger than any found in the SPS points displayed in Table [H We 
note that for each model A^^ < 10 across the mo — M1/2 plane, which suggests that these 
models do not require a precise tuning of hf. Furthermore, for models with rj^ = 1/3, we 
note that A^^ < 5 within most of the region with A < 5. 



4.2 Fine-tuning at large tan^^ 

For all of the parameter scans presented here we have chosen a fixed value of tan/3 = 10. 
For larger values of tan /3 = 20, 30 and 40 we have found that the essential structure remains 
the same; the low scale focus points still permit regions in the mo — M1/2 plane that satisfy 
the Higgs mass bound with low fine-tuning. The only significant effect of increasing tan P 
is the shift of the focus points towards lower energy scales, which is a consequence of the 
slower running of the soft mass due to the larger bottom and tau Yukawa couplings. 
As a result, is driven to less negative values, which has two implications. First, the 
magnitude of \n\ is reduced which, according to Eq. fl5U]) . also reduces |A^|. Second, the 
region of parameter space that permits EWSB is reduced, and vanishes completely for 
tan /3 ^ 50. Hence the fine-tuning tends to decrease for larger values of tan/3, however this 
comes at the cost that EWSB is more difficult to achieve. 



5 Phenomenology 



One of the main characteristics of a model with a low scale focus point are the ratios, 771 
and ?73, of the gaugino masses. From the running of the gaugino masses it is well known 
that: 



Mi(Q) 



O.Sr^i 



M.(Q) «.(g) ^'^^\ 
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where the final equahties hold for the choice Q ~ 1 TeV. Thus one can, quite directly, 
reconstruct the ratios rji and 773 by measuring the mass of these states. By comparing these 
ratios with the 'preferred region' of Fig[T], one can test the hypothesis that a focus point is 
responsible for the little hierarchy between the scales of electroweak and supersymmetry 
breaking. 

Another characteristic, which we identified in the previous section, is the existence of 
a light Higgsino with < 200 GeV. In regions of parameter space with \Mi\, M2 » \fi\, 
the lightest supersymmetric states will be a chargino and two neutralinos, which will have 
near degenerate masses. The splitting between the mass of these states is given at tree 
level by: 

m,o -m,o = Ml(± + ±]+ (39) 



^2 \Mi M2J 

- = WJ^ + ^]+lMl(^-^]e sin 2/3 + O(^) (40) 



2 ^V^i M2J 2 ^V^i M2J 'M| 
where e = Hence, at tree level, m^o < m^± < m^o. For M2, |Mi| > 1 TeV, 

this mass splitting m^± — m^o < 4 GeV. One loop corrections to these masses, which are 
dominated by third generation squark/quark couplings, generate further contributions that 
can significantly affect these splittings [35]. As can be seen by the (green) light-shaded 
region of Fig [5l for /i > there are regions of parameter space in the 54, O-II and 770 
where, for sufficiently large M1/2, the chargino is driven lighter than the neutralino as a 
result of these radiative corrections. For /i < the chargino is instead driven to a higher 
mass, and becomes heavier than the X2- The remainder of the SUSY spectrum, comprising 
the remaining gauginos and scalars, can be significantly heavier than these states, with 
masses > 1 TeV. 

If \Mi\ < \fi\, then an additional light neutralino will also be present in the spectrum. 
This occurs in the 54 and 210, for which |?7i| < 1. In these models there are regions of 
parameter space where the LSP is either dominantly Bino or Higgsino, or a mixture of the 
two. This is significant for the dark matter abundance, as we will discuss in more detail 
below. For the 210, with rji = —1/15 the Bino is very light; indeed one can estimate from 
Eq ([38]) that Mi < 45 GeV for M1/2 < 1.5 TeV. When applying the LEP limits on direct 
particle searches we placed the following requirement on the LSP: 

m^o > 46 GeV (41) 

in accordance with the PDG [22]. This limit is violated over much of the 210 parameter 



space, which is indicated by the dark shading (blue) in Figs 2(b) and 3(b) However 



as pointed out in |36], the limit of Eq fl4T]) is inappropriate for a neutralino that is light 
and dominantly bino, since it has suppressed couplings to the gauge bosons. Therefore 
it is possible that a significantly larger region of the ttlq — M1/2 is consistent with present 
observations than is indicated. This, however, requires a more detailed phenomenological 
analysis that is beyond the scope of this work. 

The presence of light Higgsino-like charginos and neutralinos has important conse- 
quences for a number of observables, in particular for b — )■ S7, and the dark matter 
abundance. We will now consider each of these in more detail. 
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5.1 6 — > S7 

FCNCs arise at one-loop within the SM, hence new physics that hes near to the weak scale 
can generate contributions at the same order. As a result, precise measurements of these 
FCNCs can provide a strong constraint on new physics models. It is well known that the 
branching ratio Bt{B — )■ Xs'~f) is particularly constraining for models such as the MSSM, 
because of its precise experimental determination and the improving accuracy of theoretical 
calculations. Combining data collected by several experiments, the current world average 
for this branching ratio is given by [37] : 

Bi(B ^ Xs-f)E-,>i.6 Gev = (352 ± 23 ± 9) X 10^^ (42) 

for a photon energy > 1.6 GeV in the rest frame of the B meson. The two errors are, 
respectively, a combined systematic and statistical error, and a systematic error resulting 
from the lineshape functioiJ^I. The theoretical calculation performed within the SM has 
been completed at NLO, and some progress has been made towards including contributions 
at NNLO. Currently the best theoretical estimate of this branching ratio at 0{a'g) in the 
SM is [38]: 

Bt(B ^ X,7)i?f>\.6 GeV = (315 ± 23) x IQ-^ (43) 

where non-perturbative effects are the dominant source of uncertainty. Supersymmetric 
models generate additional contributions to this branching ratio. These, which are cur- 
rently known up to NLO, are dominantly generated by one-loop diagrams involving a top 
and charged Higgs or a stop and chargino. 

Using the publicly available code superlso 2.4 [39], which includes both the NLO su- 
persymmetric and NNLO contributions from the SM, we have evaluated this branching 
fraction for each of the models discussed in Section HJ The results of this calculation are 
shown as the solid black contours in FigEl for Aq = 0, tan/3 = 10 and /i > 0. The case 
/X < is shown in FigO The contours correspond to the na (for n ^ Z) deviation of the 
prediction from the central experimental value, where a = 34 x 10~^ combines all experi- 
mental and theoretical errors in quadrature. We have used the error in the SM prediction 
to estimate the theoretical error. 

Because the SM prediction is la below the experimentally measured value there is a 
preference for SUSY models that generate a small positive contribution to Br (5 — )■ Xs7). 
The contributions from H^t loops are always positive, whilst the sign of loops involving 
X^t is given by [lAt [50], which from Eq (136|) is given by —firj^. 

The x^t contribution dominates in the O-II and 210 models, as a result of the small, 
positive value of rj^ which, through RGE running, tends to reduce the mass of the stops 
(relative to the H^) and give rise to significant stop mixing. As a result there is a tension 
with the experimentally measured value of Br(i? — t- ^^7) in the /i > case, since the 
X^t loops generate a large negative contribution which pushes the theoretical prediction to 
more than 2a below the experimentally measured value over most of the mo — M1/2 plane. 
Hence there is a preference in these models for n < 0, which permits the prediction to lie 

"'^^See |37j and references therein for further details. 
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Figure 5: Contour plots of the branching ratio Br(i3 — )■ Xg^) in the mo — -^1/2 plane. This is for 
the hypersurface in parameter space with tan /? = 10, = and > 0. Solid (black) contours 
are shown for Br(i? — )• Xg"^) x 10^ = 3.52 + no", where n & 'L and the error a = 0.34 combines 
all experimental and theoretical errors in quadrature. The (red) dashed and (red) dot-dashed 
contours indicate where the Higgs mass, m/jO, is 111 GeV and 114 GeV, respectively, whilst the 
dotted contours correspond to the |A| = 5 and |A| = 10 contours from Fig[2j 



within la of experiment for large values of itlq and M1/2, which coincides with the region 
that generates a sufficiently heavy Higgs. 

The 54 and 770, however, have larger, negative values of 773, and so tend to possess 
heavier stops and small stop mixing. Because of this the H'^t and x^'^ contributions are 
comparable in these models. For the > case both contributions are positive and as a 
result there is a significant region in the mg — M1/2 plane where the prediction lies within la 
of experiment. In the /i < case there is some cancellation between the two contributions 
which leads to a prediction that is not significantly shifted from that of the SM. 
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Figure 6: These are the same plots as in FigO but for the case < 0. 



The measurement of the muon's anomalous magnetic moment has been performed to such 
precision that it is now sensitive to contributions from hght supersymmetric states. Hence, 
potentially, it provides an important constraint upon all SUSY models. The latest experi- 
mental result from E821 |¥T] : 

a;fP = 11659208.0(5.4)(3.3) x 10"^° (44) 

is somewhat larger than the current SM prediction for a^, calculated using the latest 
e+e" TTTT data [12]: 

aj^ = 11659180.5(5.1) x 10'^° (45) 

where the error is dominated by uncertainties in the hadronic contributions. The discrep- 
ancy between the two results, adding all errors in quadrature, is: 

= a;fP - aj^ = 27.5(8.1) x 10"^° (46) 
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Figure 7: Contour plots of Ja^, in the mo — Afi/2 plane. This is for the hypersurface in parameter 
space with tan f3 = 10, Aq = and fi > 0. Dashed (black, with alternating long and short dashing) 
contours are shown for 5a^ x 10^^ = 27.5 + na', where n € Z and the error a' = 8.1 combines 
all experimental and theoretical errors in quadrature. The (red) dashed and (red) dot-dashed 
contours indicate where the Higgs mass, m/jo, is 111 GeV and 114 GeV, respectively, whilst the 
dotted contours correspond to the |A| = 5 and |A| = 10 contours from Fig[2j 



which is a discrepancy of 3.4(T. There are a number of unresolved issues with this dis- 
crepancy, particularly regarding the true size of the hadronic contribution!^, which has a 
substantial impact upon the theoretical prediction. Despite this, it is interesting to con- 
sider the contribution to from supersymmetric states. We have evaluated this using a 
routine included in superlso, and the results are shown as the dashed contours in Figs [7] 
and El 

The dominant contribution to (5a^ arises from one-loop diagrams involving charginos, 

^^If the hadronic contribution is calculated using data from r-decays, instead of e^e~ data, the discrep- 
ancy is substantially reduced to 1.4cr. See [13] for further discussion. 
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Figure 8: These plots are similar to those in Fig [71 except they illustrate the case < 0. 



neutralinos and slepton states. One can see in Figs [7] and |H] that the shape of the contours 
in the low mo region differs somewhat in models with 771 < (54 and 210) compared to 
those with r^i > (770 and O-II). This is the result of chargino and neutralino contribu- 
tions having signs [iMi and —\iM\ respectively; hence a partial cancellation occurs when 
M1M2 > 0. It is clear that the regions which satisfy the Higgs bound typically generate 
contributions to that are too small to account for the discrepancy observed in fH^ . This 
is caused by the relative heaviness of slepton states in these models, which is a consequence 
of the large values of either M1/2 or mo required to satisfy the Higgs bound, combined with 
our assumption of scalar mass universality. Until the uncertainties in the result are 
resolved, one cannot say whether these models are already excluded by this result. If the 
discrepancy between the SM and experiment persists then one may consider further model 
building in order to account for this. For example, one may consider introducing mild 
hierarchies into the soft masses by splitting the first two generations from the third. This 
would permit spectra with the light sleptons necessary for a significant contribution to 
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5a^j_. However, without a concrete model that can generate such large inter-generational 
splittings we will not pursue this further here. 



5.3 VL 



CDM 



The current best-fit of the ACDM model to the combined data sets of: the five year WMAP 
observations, measurements of Type-la supernovas and baryon acoustic oscillations, has 
determined the following average for the density of cold dark matter (CDM) \AA\ H5] : 



(^cdm/^^ = 0.1131 ± 0.0034 (47) 

Despite this and other compelling evidence for its existence, the precise nature and compo- 
sition of dark matter are unknown. However, the lightest supersymmetric particle (LSP), 
which is stable under the assumed R-parity symmetry, can be generated in the early uni- 
verse and, provided it is a neutralino, is a viable candidate for one component of the 
observed cold dark matter abundance. Consistency with observations and the ACDM cos- 
mology therefore requires that the neutralino abundance, fi^, is less than ^cdm- Under the 
assumption that neutralinos are only produced through the mechanism of thermal freeze- 
out, and assuming a standard thermal history of the universe, the abundance depends 
only upon the spectrum of supersymmetric states. Hence this upper limit provides an 
important constraint on the allowed regions of parameter space. 

Schematically, Vt^ oc (cxv)^^, where (av) is a thermally averaged cross-section [50] that 
combines the cross-sections for annihilation of the Xi ^Iso its co-annihilation with 
any other supersymmetric states present at the time of freeze-out. Because the lightest 
neutralino is a superposition of several gauge eigenstates: 

Xi = NfB + N^W^ + N^^H^ + NfHd (48) 

its couplings, and therefore Q^, depend critically upon the neutralino's composition, i.e. 
the coefficients Nl that are determined by the neutralino mass matrix. 

As we have already discussed, models with a low scale focus point contain a light Hig- 
gsino. In models with It^iI > 1, as in the 0-H and 770, Mi, M2 ^ over most of parameter 
space, and so the lightest neutralino becomes a (near) pure Higgsino state. Because of its 
gauge interactions the Higgsino can annihilate directly into gauge bosons, whilst the near 
degeneracy of the lightest neutralino and chargino gives rise to co-annihilation. An analytic 
calculation performed in |52] gives the abundance of a Higgsino pure state as: 



hence Q^h"^ ^ 4 x 10^^ over the entire parameter space of these models, which is safely 
within the bound given by ( H71) . Thus without positing additional (non-thermal) mecha- 
nisms of production, the Higgsino LSP can only give a subdominant contribution to the 
total cold dark matter abundance. The prospects for the detection of this Higgsino-like 
component have been studied previously in the context of non-universal gaugino masses, 
and we refer the interested reader to the literature 
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Figure 9: Contour plots of ™ ^'^'^ ^-'-0 models. These are for the hypersurface in 

parameter space with tan /3 = 10 and Aq = 0. The narrow, yellow region, which lies between the 
regions of over- and under- abundance, has a dark matter abundance that satisfies, within 3c7, 
the constraint given in Eq (j47p . The (red) dashed and (red) dot-dashed contours indicate where 
the Higgs mass, m^o, is 111 GeV and 114 GeV, respectively. 



For the alternative case that \ri\\ < 1, as in the 54 and 210, there exist regions of 
parameter space where both Higgsino and Bino are light and hence where the lightest 
neutralino is a mixture of these two eigenstates. The couplings of the Bino are such that 
it may only annihilate through t-channel exchange of superpartners, such as sleptons or 
squarks. Hence the annihilation cross-section of a pure Bino LSP becomes suppressed as 
the masses of these scalars increases, resulting in a larger abundance. As was shown by 
Arkani-Hamed et. al. |52], for slepton masses mg^ > 111 GeV the abundance is already 
forced above the limit set by Eq fHTl) . As is demonstrated in Figs [5] and [6], the constraints 
from LEP2 and Br (6 — )■ 57) already excludes most of this 'bulk region' of parameter space 
that permits these light scalars. Thus an acceptable dark matter abundance requires that 
the cross-section is enhanced either by a significant Higgsino component of the LSP, or 
through coannihilation or resonance effects. In particular, as M1/2 increases one expects to 
pass from a Bino to Higgsino dominated LSP. During this transition one will pass through 
a so-called 'well tempered' region, for which the composition of the neutralino is such that 
{av) is tuned to give agreement between and f^cDM- In order to study this in more 
detail we have calculated f2^/t^ for the 54 and 210 using micrOmegas |5T]. The results of a 
scan over the mo — M1/2 plane are displayed in Fig[9l for the parameter choice Ao = and 
tan/3 = 10. The plot is divided into three main regions, depending upon the neutralino 
abundance. Within the narrow yellow band the abundance satisfies: 

0.1029 < Q^h^ < 0.1233 (50) 

which agrees with the best fit value given in Eq P7|) to within 3a. (Note that we have 
neglected any errors in the theoretical prediction of the dark matter abundance.) The 
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remaining orange and white regions of the plot indicate an excess or under-abundance of 
dark matter, respectively. Significantly, one may note that there is a region of the mo — 
Ml/2 plane, for both the 54 and 210, in which both = ^cdm and M^o > 111 GeV are 
satisfied, with A < 10. 

There have been several previous studies of neutralino dark matter in the presence of 
non-universal gaugino masses |l6l HTl HH HQ] . In particular the case M^^Mx) = Mi{Mx) < 
Ml/2 has previously been considered in [IH], where a well-tempered neutralino, with = 
flcDM was shown to be achieved for a suitable value of M3/M1/2. The points in the yellow 



band of Fig 9( special case, justified by a GUT, of this more general study, to which 

we refer the reader for a more detailed discussion of direct and indirect detection, together 
with collider phenomenology. 

The delicate tuning of the parameters required to achieve = ^cdm, does suggest 
that some degree of fine-tuning may be present in this particular 'well-tempered' scenario. 
In order to study this fine-tuning we use the sensitivity measure of Ellis and Olive [51] . 
which defines the tuning in fi^ due to a parameter Q/i dbS'. 

A? = . (51) 

which is analogous to the fine-tuning measure used for EWSB, given in (|6]). To permit 
comparison with previous studies of A^, performed for different patterns of soft masses, 
we will utilise the following parameter set: 

G {mo,Mi/2, Ao,tan/3} (52) 

(which differs from that used for fine-tuning in Mz by treating tan /3 as a parameter instead 
of /i) and, following King and Roberts [SB], define our total fine-tuning A^ = Max{|Ap|}. 
For the 54 we find that 5 < A^ < 15 in the region which has A < 5 and M/^o > 111 GeV 
(increasing as mo increases), whilst for the 210 A*^ «i 15 . To put this into context one 
can consider the fine-tunings observed for other points in MSSM parameter space. Ellis 
and Olive determined that the CMSSM, with universal gaugino masses, has A^ < 1 and 
A^ < 5 in the bulk and stau coannihilation regions, respectively, whilst A^ ^ 200 near 
the focus point [51] . In comparison, a study of non-universal gaugino masses performed by 
King and Roberts revealed that a well-tempered neutralino, obtained by bino-wino mixing, 
has A^ ^ 28 when no relation is assumed between gaugino masses [56], which was shown 
by Birkedal-Hansen and Nelson [SS] to fall to A^ < 3, 17 (depending on tan/3) when one 
assumes a relation of the form M2{Mx) = rMi{Mx) with r taking the fixed value r = 0.6. 
Hence, in comparison, the fine-tuning A^ present in the well-tempered regions of the 54 
and 210, whilst by no means optimal, appears to be fairly mild. But whilst the narrow 
yellow band seems to satisfy the constraint of Eq f l50|) with low fine-tunings A and A^, it 
remains impossible, without further model building, to motivate why this small portion of 
parameter space should be preferred over any other. 
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6 Summary and Conclusions 



In this paper we have demonstrated how certain relations between the gaugino masses, 
present at the scale Mx, can lead to a focus point in the soft mass m'jj^ that lies near to 
the TeV scale. At this focus point has a weak dependence upon the parameter M1/2, 
which we argued in Section [2] could allow a large gaugino mass M1/2 ~ 1 TeV with a fine- 
tuning less than 10. In Section |2] we demonstrated that these ratios can be generated by 
known mechanisms of SUSY breaking. In particular, when SUSY breaking is mediated by 
gravity, we demonstrated that GUT models, through the F-term of some Higgs superfield, 
can generate ratios that will realise a low scale focus point. We have also demonstrated 
that this is possible within a class of string model, in which the supersymmetry breaking 
occurs in the moduli fields. It is particularly interesting, we feel, that the ratios required to 
realise a focus point can be motivated within the context of GUT and string constructions. 

In Section H] we computed the fine-tuning of these models numerically, taking into 
account all one-loop corrections, and demonstrated that several of these models permit a 
range of parameters for which the Higgs mass is consistent with the LEP bound of 114 
GeV, with a fine-tuning less than 5. We would like to stress that the structure of the 
one-loop corrections is particularly important in these focus point scenarios, as it tends to 
significantly reduce the levels of fine-tuning expected from a naive tree-level analysis. 

Each of these models are characterised by their spectra of gauginos and Higgsinos. In 
particular one expects a light Higgsino with mass less than 200 GeV, whilst the remaining 
gauginos and scalars can be significantly heavier, with masses > 1 TeV. The light Higgsinos 
in these models can generate significant contributions to Br (5 — Xg'j). For the case of 
Aq = and moderate tan /3 we have shown that consistency with experiment requires /i < 
in models with 773 > 0, whilst the sign of fi is weakly constrained in models with 773 < 0. 
In all of the models considered the region with an acceptable Higgs mass give very small 
contributions to a^. Thus the discrepancy 6a^ between the measured and predicted values 
of is similar to that in the SM. If this result persists, and its statistical significance 
increases, one will need to introduce further non-universality into the scalar masses to 
account for this observation. Furthermore we have also identified models which possess a 
light Bino, thus giving rise to a LSP that is a mixture of Higgsino and Bino states. For 
some region of parameter space these models possess a so-called 'well tempered' LSP, in 
which Q^h'^ is consistent with the value of flcDMh"^ required in the ACDM model. We have 
examined the degree of tuning present in these well-tempered regions, which we have found 
to be ~ 0(10). This tuning is relatively mild compared to other scenarios in the MSSM 
which achieve the correct dark matter abundance. 

To conclude, the presence of a focus point in the soft terms of the MSSM may explain 
why the superpartners are significantly heavier than the weak scale. It is a relatively 
simple hypothesis, requiring nothing more than the structure of the MSSM up to the scale 
of Mx, that can be justified within known mechanisms of SUSY breaking. Furthermore, 
by measuring the masses of the supersymmetric states, it is a hypothesis that can be tested 
in the near future. 
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Appendix A: Conventions 

Here we define the conventions that we employ throughout this paper, which are similar to 
those used within SOFTSUSY J21j. We work within the context of the R-parity conserving 
MSSM, with the soft SUSY breaking sector given by: 

Aoft = -"^i - {Mai3X"Xj + ^AijkYijk(j)i(j)j(j)k + BfiHuHd + h.c.) (53) 

where (pi and A° are scalar and gaugino fields respectively. The fields if„ and Ha refer to 
the up and down type Higgs, respectively, and we define HuHd = H^H^ — H:^H^. The 
parameter /i is the usual Higgs mass that appears in the superpotential. Without any loss 
of generality we choose our phases such that both M2 and Bfi are real and positive. In 
order to avoid constraints on CP-violating phases, we further assume that all phases are 
either or vr. 



Appendix B: Numerical Solutions of the RGEs 

In this section we summarise some of the RG coefficients, defined in Eq ( ITSl) . which are 
used in Section O 

These coefficients are solutions to the two loop RGEs [20] • They are obtained using the 
boundary conditions for the Snowmass point SPS la, which was obtained using SOFTSUSY 
[21j . The relevant quantities that define the boundary conditions are as follows: 

Mx = 2.38 X 10^^ 
giiMx) = 92iMx) = 0.721 
gs{Mx) = 0.708 
ht{Mx) = 0.494 
hbiMx) = 0.050 

hr{Mx) = 0.069 (54) 

The relevant coefficients are given in the Tables [3] and H] below, evaluated at a renormal- 
ization scale Q = 500 GeV. Note that each row contains the contribution of each gaugino 
ratio, or combination of ratios, to a coefficient. In our conventions = 1. 
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